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A bstract . We determine Saturn's radio rotation period using measurements 
made by the Planetary Radio Astronomy experiment onboard the Voyager 
spacecraft. The sidereal period deduced. 10 hr 39 min 2U sec ♦ 7 sec, is 
within the 10 hr to 11 hr range of optical periods derived from a century of 
atmospheric spot and Doppler spectroscopy observations. The radio rotation 
period is presunably that of the planet's magnetic field. We propose a 
provisional Saturn longitude convention, and we provide equations to compute a 
longitude ephemeris and to transform between the proposed system and the (10 
hr 14 min) system used irr the Pioneer 11 /Saturn encounter. In addition, we 
evaluate the degree of Iwngitude smearing which could result over the long 
term from the merging of data sets organized in this system. Finally, no 
evidence of control of the radio emission by any of Saturn's satellites is 
found . 


Introduction 

Using Voyager Planetary Radio Astronomy (PRA) data, Kaiser et al. (1960) 
recently provided the first conclusive evidence for the emission of nonthermal 
radio waves from Saturn. From an initial survey of 40 days of data, they 
observed that the Saturnian kllometric radiation (SKR) is sporadic, strongly 
and exclusively right-hand polarized, and that it occurs primarily at 
kilometer wavelengths. They further reported a provisional value for the 
rotation period of Saturn's magnetic field based on the presence of 
statistically significant modulation of the radio events and on the assumption 
of planetary magnetic field control over the event occurrence statistics. The 
period they reported, 10 hr 39 min 54 sec (+36 sec), was within the range of 
optical periods derived from Doppler spectroscopy (Moore, 1939) and from 
measurements of spots (Hall, 1877; Reese, 1971). In the present paper we 
extend and improve upon the rotation pe’^lod analysis of Kaiser et al. in two 
important ways. First, by extending the analysis interval by a factor of 
about 7, to 267 days, we are able to determine the absolute sidereal rotation 
period to within VT sec. And, second, we propose adoption of a longitude 
system adequate to the needs of Voyager and Pioneer 1 1 investigators who wish 
to organize their data in a consistent fashion at Saturn's intrinsic rotation 
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rat*. We provide an equation which generates this longitude system. In 
addition, an equation Is provided permitting translation b,>tween the 
right-handed (10 hr 14 min) system used by the Pioneer 11 Investigators during 
the September, 1979 Saturn encounter and the left-handed (10 hr 39 min 24 sec) 
system derived here. 


Observations 

The present study makes use ( r Voyager PRA data extending from the date of 
first detection of SKR, January 1, I960, throrgh September 23, I960, a 267-day 
Interval corresponding to 600 rotations of Saturn. During this period we 
observed 520 SKR events totaling 300 hr of activity. (For the purpose of this 
report we define an event as the continuous detection of SKR, lasting anywhere 
from several minutes to several hours.) Taking Into account the total 
duration of Interference-free monitoring, this corresponds to an average 
occurrence probability of 71, or an average of almost 9 events (or 30 min of 
activity) per rotation. This occurrence rate was by no means uniform over the 
analysis Interval, however, reflecting the fact that Voyager 1 (VI) was 
approaching the planet, thus Increasing the likelihood of detection. For 
example, the average occurrence probability was only 3t during the first half 
of the analysis Interval but 12t during the second half. 

We have continued to use the methods described by Kaiser et al . (I960) in 
identifying SKR in PRA frequency-time spectrograms. Briefly, where 
simultaneous coverage existed. Identifications were confirmed through 
comparison of Voyager 2 (V2) records with those from VI. That is, SKR was 
identified by virtue of Its greater Intensity and earlier occurrence in the VI 
records owing to the proximity of VI to Saturn. Where dual spacecraft 
coverage did not exist. Identifications were based on the polarization and 
dynamic spectral appearance alone. However, fewer than 1511 cf the events were 
Identified in this way. Although V2 data were used in this corroborative 
fashion, VI data were used exclusively In forming the time series from which 
the rotation period was derived because the VI data constituted by far the 
larger of the two data sets. 

A sample record of multiple SKR events on August 14, I960 Is shown in Fig. 
1. By this point in time, SKR clearly dominates PRA dynamic spectra observed 
by VI. The Jovian emission, appearing primarily in the high-frequency 
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(bottom) half of the apeotrum at 0300*0330 spacecraft event time (SCET), is 
notlcably weaker than the SKR. Unusually strong and well-defined frequency 
drifting SKR events can be seen from 0230-0330 and from 1230-1U00 SCET 
centered at about 250 and UOO kHz, respectively. In addition, a broad region 
of amorphous *^KR Is seen from 0400-0900 and from 1630-1730 SCET over the same 
frequency band. 


Procedure 

We sunmerlze In this section the steps used to derive Saturn's rotation 
period and calculate an associated uncertainty. As observed by Voyager 1, all 
events were recorded In spacecraft event time In a frame which Is moving 
relative to Saturn. In order to compute .he rotation period In the desired 
sidereal frame, the start and stop time of each event was corrected for light 
travel time between Saturn and VI. Over the duration of the analysis 
interval (267 days) the magnitude of the light time correction varied from 
23.3 to 3.8 roln. Cor/ectlon for the motion of the spacecraft In azimuth 
relative to Saturn's ascending node (motion In right ascension) was not 
considered necessary since the SKR beam Is believed to be fixed relative to 
the Saturn-sun line rather than rotating with the planet (Warwick et al., 
1981). In any event, aljustment of the event times for right ascension motion 
leads to a change In the sidereal period which Is much less than the 
uncertainty in the period. 

Following adjustment of the start and stop times for light-travel time, the 
PRA data were represented by a time series consisting entirely of 1's and O's. 

The I'a correspond to the times of SKR emission and the O's to times when 
observations were made but no SKR was recorded. The time series was analyzed 
using a method of power spectral analysis applicable to data, such as ours, 
which are sampled at Irregular Intervals (Deeming, 1 975) . A least-squares 
r..ethod was then used to refine the period found by spectral analysis of the 
time series. The data were linear least-squares fit to an equation of the 
form s a ♦ br^^, where the rotation phase, of each event was expressed 
as a function of Its Integral rotation number, r^. Rotation numbering began 
on Jan 1.0, 1980. As a first step, rotation phase and number were computed 
using the estimated sidereal period derived from the spectral analysis. Any 
residual trend in the data, that is, any tendency for the event longitudes to 
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incraas* or deoraaaa with time la then indloatad by a finite positive or 
negative value of the the least-aquarea slope, b. The best estimate of the 
sidereal rotation period is then obtained through iterative regression until b 
spproaohes zero or is within one standard deviation of b. 

Finally, an estimate of the uncertainty in the derived period is obtained 
by performing the regression many times using a range of starting values 
around that which gives the best fit. Owing to the soetter of the data in 
rotation phase, some of these analyses will yield first iteration slopes which 
are statistically insignificant (b<<>^), indicating a range of rotation periods 
which are compatible with the accepted value, fhe difference between the 
largest and smallest uuch rotation period is taken to be the uncertainty in 
the pvTiod. The uncertainty so obtained is larger than that which would be 
derived from using alone (the formal error) but it is undoubtedly a more 
meaningfull estimate. 


Results 

The result of applying the above procedure up to and including calculation 

of iiie power spectrum is shown in Fig. 2. The spectrum is constructed of 

events which occurred at a single frequency, in this case 17^ kHz, and which 

-2 -1 

exceeded a flux density threshold of 10 W m Hz normalized to a constant 
1 AU distance between the observer and Saturn. Events above this particular 
threshold were easily detected over the entire analysis interval. The lower 
panel shows a low-resolution survey extending over a wide spectral window from 
periods of 10 hr to 1000 hr. The upper panel shows a high-resolution 
spectrum, of the same data, centered approximately on the principal peak in 
the survey spectrum. The average time resolution of this spectrum is near the 
practical limit of about 65 sec. Both are plotted in terms of the standard 
deviation (<>) of the spectral power above the 'noise' level inherent in <^ach 
spectrum. 

The survey spectrum (bottom panel) spans the revolution periods of all of 
Saturn's major satellites and the expected period, as viewed from Saturn, of 
persistent solar wind features such as Interplanetary magnetic field 
reversals. In greater detail, the high resolution spectrum (top panel) spans 
the ''' 10 hr to ''' 1 1 hr range of periods of Saturn's surface features and 
clouds (Moore, 1939; Reese, 1971). Examination of the survey spectrum 
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indicates only one significant peak, near 10 hr, which is shown by the 
high-resolution spectrum to occur at 10.66 hr. This feature is a factor of 
5 above the second strongest peak (near 260 hr) in the spectrum and 
approximately 16<> above the survey spectrum noise level. We attribute this 
spectral feature to the modulation of the SKR by Saturn's rotation. With 
regard to ix>ssible satellite modulations, we find no evidence of spectral 
power at the appropriate satellite periods or their heterodynes with the 
rotation period. In fact, none of tt' secondary peaks exceed above the 
survey spectrum noise level and, further, their power levels have declined 
steadily with the accumulation of data. Thus at this stage cf the analysis we 
have found no Saturn analog of the well-known phenomenon whereby some of 
Jupiter's sporadic radio emission is controlled by lo. There is also no 
indication of spectral power at the low-frequency end of the survey spectrum 
where longer-te^m modulations such as those due to solar influences might 
appear. However, far too few cycles of such long-term modulations are 
included in the time series to yield a definitive test. For example, only 
about 10 cycles of a 26-day solar periodicity could have occurred over the 
analysis interval. Finally, although not shown, no significant spectral power 
was found at periods <10 hr except at the harmonics of the rotation period. 

Looking at the high-resolution spectrum, we note the 5*’ spike cn the 
low-frequency shoulder of the SKR spectral peak. This spike occurs at a 
period of about 10.69 hr and appears to be due to a low intensity component of 
the SKR which reappears at a different, slightly longer period than that of 
the principal source. This component Is not expected to significpntly bias 
the principal spectral peak because it lies well outside the 65-sec half-power 
spectral bandwidth of the peak. This component of drifting SKR is discussed 
again below where it appears explicitly in Fig. 3> 

We have refined this power spectrum determination of the period according 
to the linear least-squares method described previously and illustrated now in 
Fig. 3. Here the SKR events are indicated by vertical bars plotted as a 
function of longitude and rotation number. (Longitude is defined arbitrarily 
in this figure in order to place moat of the the events near 180° for 
clarity.) As in Fig. 2 we have plotted only events at 17^ kHz which exceeded 

«, 1 Q -2 -1 

10 W m Hz . It is evident that one range of longitudes is favored, 
namely that from about 100 to 270 degrees; this is particularly apparent in 
the occurrence probability histogram plotted on the right of Fig. 3. Fig. 3 
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(bottom) actually Ulustratas the final organization of the data in rotation 
phase following several least-squares iterations. That is, events were first 
org>^nized in rotation phase and rotation number using the provisional power 
spectruni period (10.66 hr). This yielded a small but statistically 
significant negative value for the least-squares coefficient, b (the slope), 
thus indicating a general trend in the event longitudes toward decreasing 
values with time. Clearly the optimum period is that one which removes any 
linear drift in the event longitudes. This period was found through 
successive iteration to yield the results shown in the Fig. Using a 
sidereal rotation period of 10 hr 39 min 2U sec (corresponding to a rotation 
rate of 810.76 deg/day), the least-squares fit, shown via the dashed line, 
gave zero slope within ^ Least-squares analysis using initial starting 

periods ranging from 10.65<<5 through 10.6585 hr yielded statistically 
insignificant trends in longitude from which, as descibed in the previous 
section, we derive an uncertainty of ♦dO. 6585-10. 65^5)/2, or ♦? sec, in the 
measured period. 

Along the top of Fig. 3 is shown the time history of the rotation period as 
a function of the end point of the analysis Interval. Rotation periods were 
computed in 30-rotation Increments with each successive computation including 
all of the preceding data. Because the individual determinations are not 
independent they give a deceptively good impression as to the magnitude of the 
scatter with time; however, the figure does show how the measurements are 
converging with the accumulation of data. 

It is probably apparent from the dashed-line Intercept in the figure that 
in performing the least-squares computation only the events within the main 
peak were used, that is, those events which occurred between 100 and 270 
degrees longitude. This had the effect of eliminating 27f of the data from 
the final step of the analysis. The reason for this restriction is that we 
wished to confine the analysis to those events which presumably resulted from 
successive passages of one and the same emission cone past the spacecraft. In 
other words, rather than allow the least-squares determination to be 
influenced by radio events which may have occurred in association with side or 
back lobes cf the main beam, we restricted attention to the main-beam events 
alone. Although the problem is conceptually simpler whjn expressed and 
performed in this way, it has the disadvantage of requiring somewhat 
subjective decisions regarding the choice of main-peak longitude boundaries. 
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Thus, In order to jssure the precision of the period detcrnlnstlon, we hsve 
tested the sensitivity of t’>ie technique to s brosd rsnge of possible 
vcrlstlons In methodology, Including the use of s rsnge of longitude 
boundaries, the use of Intensity weigh :ed as opposed to unweighted 
regressions, the exclusive use of event start tltaes Instead of both start and 
stop tines, and the use of events at numeroua other frequencies In the SKR 
band. Except for an analysis at 385 kHz which differed from the proposed 
period by 1.6 none of the variations yielded s rotation period which 
differed by more than 0.7 from that quoted. Further, we have subjected the 
entire analysis procedure to an Independent test by measuring Jupiter’s 
sidereal rotation period. Because Jupiter's period Is a known quantity, this 
constitutes a test of the accuracy of the procedure. We used Planetary Radio 
Astronomy detections of Jupiter near 1 MHz made from Voyager 1 over the period 
1977 November to 1978 June (Kaiser et al., 1979), approximately equivalent to 
the Saturn data set both In terms of the extent of the analysis Interval and 
the total number of events. Using the identical prccedure followed for 
Saturn, except of course for the light time and right ascension corrections 
which applied to the Voyager-Juplter instead of the Voyager-Saturn geometry, 
we derived a rotation period of 9 hr 55 min 29.6 sec ♦ 1.1 sec, or only 0.1 
sec less than the accepted period (Riddle and Warwick, 1976). The j^l.1 sec 
uncertainty is smaller In this case than that quoted above for Saturn because, 
as noted by Kaiser et al . (1980), the Jovian emission Is more highly organized 
In rotation longitude than Is SKR. 

We noted with regard to Fig. 2 the appearance of statistically significant 
spectral power at 10.69 hr, corresponding to a second, weak component of SKR. 
This component is visible in Fig. 3 as those SKR events which are clearly not 
locked In rotation phase but which 'drift through' the plot from early to late 
longitudes. Note especially rotation numbers 225 through 300 and 350 through 
425. These events are much more conspicuous In plots in which no Intensity 
threshold has been set, implying that they correspond to a weaker SKR 
component. They art perhaps analogous to the Jovian nKOM emission which 
Kaiser and Desch (1980) have shown is related to an emission source which 
falls to rigidly corotate with the Jovian magnetosphere. This component is 
being Investigated further as to other possible distinguishing characteristics. 
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Longitude Convention 


To fsollltate coneistent use of the rotation period, we propose adoption 
of the following Saturn longitude ayatain which will be referred to as SLS 
(Epoch, 1960.0) or, simply, SLS. Like Jupiter's System III, SLS is 
left-handed (positivo west), that is, the longitude increases with time as 
viewed by a stationary observer. Also like System III, the central meridian 
is referred to the geoaietrioal rather than to the illuminated disk of the 
planet so that solar phase angle ia irrelevant. The ayatem epoch is 0 hr U.T. 
on Jmnuary 1.0, I960 at which time the prime meridian (0° sls) is defined to 
be coincident with Saturn's vernal equinox, or the ascending node of Saturn's 
orbit on its equator (Sturms, 1971). (We have chosen a system epoch of 
1960.0, rather than the more conventional 1950.0, so as to minimize any 
longitude shift resulting from future updates to the rotation period. Note 
that a change in the period of as little as 7 sec (the present uncertainty) 
would introduce a longitude shift of over 1600° over the 30-yr span between 
1950 and I960.) No correction is necessary for equinox precession, owing to 
its negligible magnitude relative to the present (and anticipated) uncertainty 
in the measured rotation period. For the same reason the small difference 
between U.T. and ephemeris time is entirely negligible. 

Lo.igltude computation is as follows. If t is the time of any given 
observation, t^ the system epoch defined above, D the observer-Saturn distance 
in km, and R the observer right ascension (angle in degrees in Saturn's 
equatorial plane) from Saturn's vernal equinox to the observer, then the 
longitude, of the meridian containing the observer is given by 

^SLS * (t-tjj - 3.86x10"’’ D) - R, (1) 

where the longitude is that facing the observer at the time the light signal 
left the planet (the central meridian longitude). Note that because the 
rotation rate constant, 810.76, is in degrees/day, t-t^^ must be expressed in 
days and fractions thereof. 

The system proposed here is directed principally at the use of Voyager 
data; however, for those who wish to order Pioneer 11 /Saturn encounter data in 
the present system, we provide the following transformation equation which is 
expressed in terms of the sub-Pioneer longitudes based on the 10 hr 14 min 
rotation period. The relationship is 
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'SLS • \ - ’*397 At. 
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where is the S1.S (1960,0) sub>Ploneer 11 longitude, !■ thn sub-Pioneer 
11 longitude using the 10 hr 14 min rotstion period (Sturms, 1971), end ^t is 
the elspsed time in hours from 12.0 hr spsoeorsft event time on 31 August 
1979. Note thst the 10 hr 14 min system Miopted for the Pioneer/Saturn 
encounter is a right-hand system, that is, the longitude decreases with time. 

Longitude Smearing 

Unlike the Jovian case in which the period is known to within an accuracy 
of at least 0.02 sec, the 7-sec uncertainty in Saturn's period may yield 
measurable smearing, or 'smoothing', of data organized in rotation phase if 
the data are collected over a significant span of time. The amount of 
smearing is straightforward to evaluate, and we do so here in the context of 
the Voyager-Saturn encounter jince a principal use of thi.i period and 
longitude system will be by Voyager-spaceoreft investigators. Propagation of 
the *7 S'io uncertainty in the period between the system epoch on Jan. 1.0, 

I960 and the Voyager 2 encounter with Saturn on Aug. 27, 1961 yields no more 
than 91° of rotation phase (longitude) error. Between the Voyager 1 and 2 
encounters with Saturn (Voyager 1 encounter is on Nov. 12, 1960), 40° of 
rotation error can accixnulate. And, over the approximate 6-day bow 
shook-to-bow shock durations of the encounter* periods themselves, no more than 
0.9° of longitude smearing can take place. This is to be compared with the 
2^0° error which would result over the same 6-day interval if the 10 hr 14 min 
period adopted for the Pioneer 11 encounter with Saturn were used. Finally, 
prope:gation of the sec uncertainty in SLS from the System epoch back to the 
Sept. 1, 1979 epoch of the Pioneer 1 1 /Saturn encounter yields a maximun error 
of 18° in calculations of using eq. (2). 

Summary 

A power spectrum and linear regression techniq>>r has yielded a value for 
Saturn's sidereal rotation rerio^ of 10 hr 39 min 24 ♦ 7 sec. This 
corresponds to a rotation rate of 810.76 * 0.15 deg/day, which we take to be 
the rotation rate of Saturn's magnetic field. A left-handed (positive west) 
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Saturn longltuda ayatam (SLS, 1960.0) la daaorlbad In aquation Tom. Tha 
convantion haa Ita prlma marldian (0° SLS) oolnoldant with Saturn* a varnal 
aquinox at 0 hr u 7. on Jan. 1, I960. Aaida from tha powar appaaring at 
Saturn* a rotation parlod and a minor ooreponant at a alightly longar pariod, 
tha wida-window powar apaotrun ia faaturalaaa indioating that thara ia no 
avidanoa for modulation of 8KR by any of Saturn*a aatallitaa. 
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Figure Captions 


Fig. 1. Voyager 1 frequenoy-time dynemlo spectrum showing Saturn kilometrio 
radiation (SKR), Jovian emission and solar type 111 bursts. Increasing 
darkness is proportional to increasing signal strength. The spectrum covers 
the 24->hr period beginning at 0 hr spacecraft event time (SCET) on August 14, 
I960. The SKR frequency extent is from about 80 kHz to almost 750 kHz. At 
this time VI was about 0.8 AU from Saturn and about 3>4 AM from Jupiter. 

Fig. 2. Low-resolution (bottom) and high-resolution (top) power spectra of 
SKR time series. Only events at a single frequency (174 kHz) and above a 
constant flux threshold were used. Each spectrum shows only a single major 
peak, corresponding to modulation of SKR at Saturn's rotation period. 

Fig. 3. A plot (bottom panel) of the same SKR data as in Fig. 2 in terms of 
an arbitrary longitude system and rotation number. The (iata are shewn after 
final organization at the adopted sidereal period as indicated by the 
least-squares fit (dashed line). The figure in the top panel shows how the 
period has varied with the accumulation of data. A scale relating calendar 
date to rotation number is shown between the two panels. The longitude 
histogram on right shows how the data are organized in rotation phase. 
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